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When a source emits a gravity-wave (GW) pulse over a short period of time, the leading edge of 
the GW signal is redshifted more than the inner boundary of the pulse. The GW pulse is distorted 
by the gravitational effect of the self-energy residing in between these shells. We illustrate this 
distortion for GW pulses from the final plunge of BH binaries, leading to the evolution of the GW 
profile as a function of the radial distance from the source. The distortion depends on the total GW 
energy released e and the duration of the emission r, scaled by the total binary mass M. The effect 
should be relevant in finite box simulations where the waveforms are extracted within a radius of 
< 10 2 M. For characteristic emission parameters at the final plunge between binary BHs of arbitrary 
spins, this effect could distort the simulated GW templates for LIGO and LISA by a fraction of 10 -3 . 
Accounting for the wave distortion would significantly decrease the waveform extraction errors in 
numerical simulations. 



PACS numbers: 

I. INTRODUCTION 

The observation of gravitational waves (GWs) is ex- 
pected to open a new window on the universe within the 
following decade. First generation GW detectors (In- 
LIGO |zg, VIRGO [z3, TAMA |z|, GEO [zjj) are al- 
ready operating at or close to their design sensitivity lev- 
els and the development of the next advanced-sensitivity 
GW detectors (Advanced LIGO [HJ, Advanced Virgo 
[HI, LCGT [13) and the space-detector LISA H are 
well underway. It is now increasingly important to fully 
understand the precise characteristics of the GW wave- 
forms that we expect to observe. 

The most luminous GW sources are expected to be 
associated with mergers of BH binaries. The physical 
understanding of these sources has greatly improved by 
recent breakthroughs in numerical relativity [l], [U, [I[ . It 
is now finally possible to simulate the merger of a BH bi- 
nary, from the initial circular inspiral, through the plunge 
to a common surrounding horizon, to the final ringdown, 
as the remnant settles down to a quiescent stationary 
Kerr-BH. It is widely believed now that existing simu- 
lations are sufficiently precise to allow targeted searches 
for these waveforms in real data In fact, it has been 
recently shown that the errors are not even limited by 
the numerical precision of the simulation (~ 10 -5 ), but 
the GW extraction method itself entails a much larger 
uncertainty (~ 10~ 3 ) [5|. In this paper, we demonstrate 
that the self-gravity during the propagation of gravita- 
tional radiation in the zone of wave extraction of numer- 
ical simuations leads to the distortion of the waves, cor- 
responding to similar magnitude modifications in typical 
cases. 
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A. Description of the effect 

Let us imagine a compact spherically-symmetric con- 
figuration of matter (representing the remnant) and a 
rapidly expanding sphere of massless particles (represent- 
ing the radiation) carrying away some of the initial mass 
of the system (Fig. Q]). First, let us assume Newtonian 
gravity and spherical symmetry. In this case, the var- 
ious shells are pulled back only by the gravity of the 
mass interior as if it was concentrated to a point mass 
at the center of the sphere, and the effect of the outer 
enclosing shells exactly cancels out. Thus, the particles 
on the outermost shell are always attracted by the total 
mass, including the mass of the radiation, but the inner- 
most shells experience only the gravity of the remnant. 
Therefore, the gravity of the radiation implies that the 
innermost shells of radiation will be continuously catch- 
ing up to the outer boundary during their journey from 
the source to the observer. 

Do we expect an analogous effect to exist also for grav- 
itational waves in full general relativity? First, let us 
consider conventional (i.e. non-gravitational) radiation. 
In analogy to the Newtonian gravitational pull, relativis- 
tic test particles are slowed down by gravity: the null- 
geodesics in a gravitational field experience the so-called 
Shapiro delay @, decreasing the radial coordinate ve- 
locity with increasing gravity. Furthermore, according 
to Birkhoff's theorem, the spacetime outside a spheri- 
cally symmetric distribution of energy is equivalent to 
the spacetime of a point-mass placed at the center of 
the sphere, the Schwarzschild metric, and the spacetime 
inside a cavity is the free-space Minkowski spacetime. 
More generally, the spherically symmetric expansion of 
collisionless radiation is a known simple exact solution 
of the Einstein equations, the Vaidya metric @,H[ . This 
solution has exactly the same characteristics as the New- 
tonian example, whereby various shells react only to the 
mass interior to them, i.e. they move on world lines ne- 
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FIG. 1: A sketch of the effect under consideration. The co- 
alescence of two BHs in a binary of total initial mass Mq 
results in the emission of a burst of gravitational radiation 
which carries away a non- negligible e fraction of Mo- The 
remnant BH mass is Mf. The proper temporal width of the 
wave-packet for a hypothetical observer fixed at a radial dis- 
tance r is At. As the packet propagates outwards, it (1) 
expands due to gravitational redshift of the initial mass Mq 
(solid lines) , (2) contracts due to the mean self-gravity of the 
radiation (dotted line, At'), and also (3) distorts its profile 
due to the self-gravity of the radiation (not shown) . As a re- 
sult the inner shells begin to catch up, and the proper time 
separation in excess of gravitational redshift from the front 
of the burst decreases with distance. Consequently, the net 
luminosity of the radiation burst changes with distance. 

glecting the exterior shells and the effect of the interior 
shells is the same as if they were concentrated to a point 
mass at the center. 

Next, let us turn to the case of gravitational radiation. 
The effect of the self-energy of gravitational radiation can 
be accounted for in the first nonlinear-order approxima- 
tion of the Einstein field equations by attaching terms 
of order h 2 to the stress-energy, T^ , considering these 
terms as sources in addition to the regular radiation fields 
[UBS E BO- Here hi, is the wave amplitude 
which is the correction to the background metric. If the 
wavelength of the GW wave-packet is much smaller than 
the size of the wave envelope, the evolution of the wave- 
packet is determined by the WKB cycle-averaged effec- 
tive stress-energy tensor [lfj . , independent of the spe- 
cific wave-characteristics of the radiation. In this regime, 
we may treat the GW packet as an ensemble of relativis- 
tic particles for which our previous arguments apply. In 
conclusion, we anticipate that 

(i) the wave-envelope will continuously expand due to 
the redshift of the initial mass of the binary, 

(ii) it will contract due to the self-gravity of the radia- 
tion, and 

(iii) in analogy to electrodynamics, we expect that the 
distortion of the wave envelope would lead to a 



continuous adiabatic modification in the GW fre- 
quency. 

The purpose of this paper, is to quantify these expec- 
tations for typical BH merger waveforms using simple 
models and to demonstrate that this effect should be 
accounted for in relation to numerical simulations and 
observed merger waveforms. 

B. Related literature 

To our knowledge the effect of self-gravitational dis- 
tortion of GWs had not been explicitly recognized previ- 
ously. We elaborate on the relation of the self-distortion 
effect to numerical general relativity, analytical investiga- 
tions like the multipolar post-Minkowskian (MPM) and 
post-Newtonian (PN) theory, and the studies of the scat- 
tering of gravitational radiation in curved spacetimes. 

The self-gravitational distortion of GWs is a relatively 
small effect on short scales currently accessible to numer- 
ical simulations. Current state-of-the-art simulations of 
binary BH mergers are restricted to the central strong- 
gravity domain near the black holes, and extract gravity 
waves from the boundary of this domain. The standard 
method of extracting and extrapolating the waveforms to 
larger distances, is based on the Regge-Wheeler-Zerilli- 
Moncrief perturbation formalism [lj| US, (3 ■ This is a 
linear-order representation of the Einstein field equations 
and so it neglects self-energy effects of order h 2 . Cumula- 
tive nonlinear effects like the self-distortion effect should 
lead to systematical deformations of the linear waveform 
extracted at different radii, which can in principle be dis- 
covered by a rigorous convergence test. In fact, nearly 
all papers on simulated merger GWs study the conver- 
gence behaviour in some detail. However, due to com- 
putational limitations, the extraction radius is currently 
restricted to r < 50M, and the extraction has been pre- 
formed on only a few, typically 3-4 different radii with 
the extrapolation done empiricall y b ased on these radii 

si s, s i, m m, m m Hi m m m ^ . Recently, 

Pazos et al. Q reported a systematic effect of order 10 -3 
(for extraction radii r < 80 AT), which is much larger than 
the numerical precision of the simulation ~ 10~ 5 . This 
is roughly the same order of magnitude systematic effect 
that we expect for the GW distortion for the given radii 
(see below). Note however, that Berti et al. [25| showed 
that the convergence behavior is also largely sensitive to 
simulation assumptions. On small scales probed by these 
simulations, other near- field nonlinearities might also be 
of equal importance. 

A precise treatment of the waveforms at large radii 
is possible by analytical methods, such as the MPM- 
expansion introduced by Thorne [141 ] (see als o 11311) and 
developed further by Blanchet & Damour [H, I29ll30l| . In 
Ref. [14| . Thorne introduced the concept of "local- wave 
zone," which is outside the dynamical zone of wave gener- 
ation, but where nonlinear effects are still important. In 
this region the propagation of gravity waves is expressed 
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in terms of an expansion in powers of hij as an infinite 
sum of multipole contributions with rapidly decreasing 
amplitude, which needs to be matched to the dynamical 
gravitational field generated by the source. The GWs 
in the local wave zone are given formally by the MPM 
expansion, whose terms correspond to different powers 
of the gravitational coupling constant G. In the PN ap- 
proach, the dynamical wave generation is calculated an- 
alytically in an infinite series in the inverse speed of light 
c~ 2 . Matching the PN and MPM expansions in the local 
wave zone is a successful method for the calculation for 
steady source GWs produced by relatively slow motions, 
like the inspiral phase of BH mergers where the distance 
between the BHs is large enough to allow an adiabatic 
quasicircular orbit at r > 6M. To date, the PN wave- 
forms for circular binary inspirals are available to 3.5PN 
order for general mass ratios [3l| (which is the highest 
order that is expected to have a measurable contribution 
for circular inspirals by a LIS A- type detector (32J) and 
5.5PN order for extreme mass ratios and no BH spins [33| . 
To our best knownledge, the self-gravitational distortion 
effect has not been identified in these works. However, 
in this paper we show that the radius for a fixed lumi- 
nosity shift or a fixed frequency shift is linearly sensitive 
to the energy density of the radiation. The luminosity at 
the inspiral phase of binaries is less than 1% of the lu- 
minosity at the final plunge 26]. Therefore, even if it is 
negligible for inspirals, the self-gravitational modulation 
of GWs could be significant for the final plunge. 

We expect the self-gravitational distortion of GWs to 
be consistent with the PN/MPM expansion, and to have 
corresponding PN tail counterparts 34]. The tails of 
GWs are caused by the scattering of linear waves on 
the spacetime curvature generated by the total mass- 
energy of the source [35| , which is related to the Shapiro 
time-delay [6j of the radiation crawling out of the back- 
ground gravity of the source. At 2.5PN order beyond the 
Newtonian quadrupole formula, GW tails scatter off the 
monopole field of the remnant [36], [37., 38] . Furthermore, 
above 3PN order, the tails of the tails are produced by 
curvature scattering of the tails of the waves themselves, 
associated with the cubic nonlinear interaction between 
two mass monopole moments and the mass quadrupole of 
the source [39| . In this paper we show that the modifica- 
tion caused by the spherical self-gravitational distortion 
of GWs is to lowest order proportional to the original 
waveform (i.e. without this effect) times the energy den- 
sity of the waveform. Since the energy density is pro- 
portional to the square of the amplitude, this possibly 
implies to lowest order a monopole-quadrupole 2 type in- 
teraction counterpart. 

The self-gravitational distortion effect is also related to 
the "memory effects" (or "hereditary effects") of gravi- 
tational radiation, since it is a cumulative effect that de- 
pends on the full past history of the radiation, as opposed 
to regular PN terms which only depend on the instanta- 
neous retarded fields. Other known hereditary GW ef- 
fects are the tails of GWs [30l [40| . the Christodoulou 



effect jUSi], and the GW recoil kick [H|. 

Finally, nonlinear effects were also examined for the in- 
teraction of pla ne GWs on a free space (i.e. Minkowski) 
background [H, [H, [H, H3, H El ■ The nonlinear terms 
in the scattering problem are found to exactly cancel up 
to fourth order, leaving no sclf-phase-modulation effect 
for GWs in vacuum to this order. However, the geome- 
try of this case is very different from the one discussed in 
this paper where a curved background is initially present. 
It is also unclear whether there are self-phase modulation 
effects at higher nonlinear orders. Fortunately, our ap- 
proach does not face such convergence issues, since we 
adopt the exact (i.e. non-perturbative) solution of the 
Einstein equations in the spherical WKB approximation 
of an expanding radiation shell. 

The present paper aims to quantify the self- 
gravitational distortion effect for recently compiled 
merger waveforms. In § HH we list the main properties of 
the waveforms relevant for our study. In § 11111 we present 
our analysis in spherical symmetry, and derive the results 
for the self-gravitational distortion of signal duration and 
the luminosity profile. In § IIVI we summarize the main 
conclusions, and then discuss their implications. Finally, 
we discuss the validity of the spherical approximation 
and consider the possible effect of the anisotropy in the 
Appendix. We use units with G = c = 1 and a metric 
signature of (—1, 1, 1, 1). 



II. MERGER WAVEFORMS 

To illustrate our effect, we adopt a simplified treatment 
for the merger GW waveforms. Table [T] lists the total ra- 
diated mass Arritot relative to the initial total mass Mq 
for various encounters between compact objects found in 
the literature. While all of these encounters would have a 
nonnegligible GW self-gravitational effect, the inspiral- 
merger-ringdown events are expected to have the most 
prominent event rates for interferometric GW detectors. 
Inspiral detection rate estimates are between 0.3 — 3yr _1 
and several per day for NS/BH mergers for inLIGO and 
adLIGO, respectively [H(J; 3 and lOOyr -1 for stellar 
BH/BH inspirals in globular clusters [51| and in galac- 
tic nuclei [52J with adLIGO, respectively; 1-100 yr _1 and 
30 yr -1 for supermassiye (SMBH) and intermediate mass 
BH (IMBH) [H, H, i| and SMBH/SMBH inspirals 
[56|,[57j with LISA, respectively. The dynamic time of the 
encounter is proportional to the total mass; consequently 
we do not consider extreme mass ratio inspiral-mergers 
(EMRI) , as the GW luminosity of these sources is much 
smaller. Although high velocity stellar BH encounters 
can be very bright in GWs, these events are expected to 
be rare, less than lyr -1 for adLIGO or LISA [11]. In 
this analysis, we focus on equal mass inspiral-merger- 
ringdown waveforms. 

Generally, the waveforms can be expanded in multi- 
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TABLE I: Total radiated mass relative to the initial total 
mass, e, for bright GW encounters. For a detailed comparison 
of BH inspiral computations see Ref. 
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Encounter 11 
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Kers. 
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( s± s 2 \ 
V"^7' ™?/ 


4-4-- b 

tat ion 




[%] 


NS-BH 


head-on 


(0.0,0.0) 




[60] 


0.01 


BH-BH 


head-on 


(0.0,0.0) 




[6JJ 


0.05 


BH-BH 


head-on 


(0.1,0.1) 




[61] 


0.06 


NS-BH 


orbiting 


(tidal c ,0) 




[62, 63] 


0.1 


BH-BH 


head-on 


(0.2,0.2) 




[61] 


0.12 


BH-BH 


head-on 


(0.0,0.0) 




[64] 


0.13 


BH-BH 


whirl 


(0.0,0.0) 




[24] 


0.5-3 


BH-BH 


grazing 


(0.9,0.7) 


- _L - 


f65] 


0.9 


BH-BH 


grazing 


(0.9,0.7) 


+ -L + 


[65] 


1.0 


BH-BH 


grazing 


(0.0,0.0) 




[65] 


1.2 


BH-BH 


inspiral 


(0.2,0.2) 


+ 1 + 


[20] 


1.8 


BH-BH 


inspiral 


(0.1,0.1) 


- 1 - 


[20] 


2.0 


BH-BH 


inspiral 


(0.2,0.2) 


- 1 - 


[20] 


2.0 


BH-BH 


inspiral 


(0.8,0.8) 


- 1 - 


[22] 


2.2 


BH-BH 


inspiral 


(0.1,0.1) 


+ 1 + 


[20] 


2.4 


BH-BH 


inspiral 


(0.0,0.0) 




[20] 


2.5 


BH-BH 


inspiral 


(0.0,0.0) 




[3, 21J 


3.2 


BH-BH 


inspiral 


(0.0,0.0) 




[19, 22] 


3.5 


BH-BH 


inspiral 


(0.0,0.0) 




[25, 27] 


3.7 d 


BH-BH 


inspiral 


(0.1,0.1) 


+ II + 


[26] 


5.2 


BH-BH 


inspiral 


(0.8,0.8) 


general 15 


[23] 


5-6 


BH-BH 


inspiral 


(0.8,0.8) 


+ 11 + 


[22] 


6.7 


BH-BH 


r. whirl 


(0.0,0.0) 




[24] 


15'' 


BH-BH 


r. head-on 


(0.0,0.0) 




[66] 


16 


BH-BH 


r. head-on 


(1.0,0.0) 


1 


[67] 


17 _n_g 

Ll 1/4 


BH-BH 


whirl 


(0.0,0.0) 




[24] 


24' 1 


BH-BH 


r. whirl 


(0.0,0.0) 




[24] 


100'' 



"Here "head-on" stands for a direct collision with v \\ r initially, 
"orbiting" stands for the tidal stripping of a NS by a BH in close 
orbit, "grazing" stands for inspiral collisions in which the initial 
separation is within the final orbit in a merger, "inspiral" is the 
complete inspiral-merger— ringdown event, "whirl" stands for parti- 
cles approaching from infinity with some impact parameter leading 
to a quasi-circular whirl-type orbit before merger, "r. whirl" and 
"r. head-on" corresponds to relativistic initial velocities » rj 1 at 
r > M. 

'Here we give the sign of Si ■ J, the relationship between Si and 
S2 , and the sign of S2 ■ J, where J is the orbital angular momentum. 
In case of a head-on collision with spins, we give the initial direction 
of the spin relative to the separation vector. 

C NS tidally locked 

d Ref |25H provides the results for different mass ratios, mi/ni2 = 
1—4, and found that Am <x r) 2 , where r) = mim^Kmi + m^) 1 < 
1/4. 

E 8 different choices of spin orientations 

■^In case the impact parameter is small enough to end up in a 
merger. 

9 Calculated for mi -C mi- (See r) above at d .) 
h ln case the impact parameters are fine tuned for the binary to 
approach the unstable circular orbit. 
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where (f> is the high frequency GW phase which is related 
to the instantaneous frequency through / = d0/di; A^ u 
is a slowly varying envelope describing how the instanta- 
neous amplitude changes over the waveform, the index n 
labels the various polarizations and multipoles, and d,L is 
the luminosity distance, and z is the redshift. For binary 
inspiral merger waveforms, the (I = 2, m = 2) multipole 
(i.e. quadrupole) dominates the waveform [25L 125. Y2^\ . 
For the sake of simplicity, we restrict our attention to a 
single monochromatic wave. 

During a single GW cycle, corresponding to a time 
interval / , the envelope can be regarded as con- 
stant. In the WKB approximation, the energy carried by 
the radiation can be calculated as a cycle-averaged quan- 
tity. The effective stress-energy tensor of the radiation is 



rpfj.v 



1 A 2 



k"k u , 



(2) 



where A 2 — A^A^ is the squared effective amplitude 
and k^ — (j)^ is the wave number. 

Each infinitesimal volume of the wave can be at- 
tributed the mass-energy that it carries according to Eq. 
(|2j). In the spherically symmetric approximation, the to- 
tal luminosity (or more precisely the graviton number) is 

El 



L = 



dE 
drdt 



^r 2 T l 



A 2 
16' 
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With this equation it is possible to obtain the waveform 
[A(t),f(t)] within radius r from the luminosity function 
L(t). 

Based on the waveforms derived by numerical simula- 
tions (such as Fig. 25 in Ref. :26[), we adopt the following 
simple fit to the effective luminosity 



L(t) = L 




if - t < t < 
if < t < ti (4) 
if t > ti 



where the intervals i<0,0<<<£i and t\ > t corre- 
spond to the late inspiral/final orbits, the peak luminos- 
ity at the plunge, and the ringdown phases, respectively, 
—to sets the initial time of the calculated profile, t\ rep- 
resents the characteristic timescale of the most intensive 
part of the radiation, ti sets the ringdown decay rate, and 
Lo is the normalization luminosity. Numerical simula- 
tions show [1^1 that the characteristic frequency of the ra- 
diation rapidly increases after the inspiral and saturates 
at bj = 2nf ~ A1q 1 /2 at t ~ 0, where Mo is the initial 
mass of the source. The characteristic number of wave 
cycles during the brightest phase is N — At/f^ 1 ~ 3 
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over a time At = 12ttMq. Note that we assume that 
the WKB method is applicable to the waveform, imply- 
ing that Lit) does not change greatly over a single cycle. 
This condition is just marginally satisfied for these wave- 
forms. 

For a simple analysis, we assume that the total lu- 
minosity crossing a sphere at infinity L^ft) is given 
by Eq. (fj| with the following parameters: (io>ii>*2) = 
(100, 10, 5)Mo, and use At = 12ttMq in reference to av- 
erage quantities below. Under the WKB approximation 
the luminosity (j4|) evolves independently of the carrier 
frequency / ~ 1/ [AttMq). The L$ normalization is set 
using the total radiated mass of the system as a frac- 
tion of Mo by e = Amtot /Mq for the inspiral-mergers ([3|) 
listed in Table HI For pedagogical comparison purposes 
we distinguish the inspiral events only based on the nor- 
malization Lq, and do not consider the variations in the 
shape of the waveform (e.g. At). 



III. QUANTITATIVE ESTIMATES IN 
SPHERICAL SYMMETRY 



To give a first quantitative estimate of the magnitude 
of the wave-packet distortion due to self-energy, we start 
by computing the propagation of unpolarized radiation 
packets in the spherically symmetric Vaidya spacetime. 
The possible effect of anisotropics is discussed in the Ap- 
pendix. 

The Vaidya metric [7J in radiation coordinates 
(u, r, 9, <p) is 



ds 2 



1 



2m(u) 



du 2 



2drdu 



r 2 dfl 



(5) 



de 2 



This metric is an exact 



where dfl 

(i.e. nonperturbative) solution of Einstein's equations in 
spherical symmetry in the eikonal approximation to a ra- 
dial flow of unpolarized radiation. Here u is the retarded 
time parameter which is constant along the world lines 
of radially outgoing radiation, and m(u) describes the 
mass function interior to u. Outside the radiation (i.e. 
where m(u) is constant in the space-time), the Vaidya 
metric ([5]) is the Schwarzschild solution in Eddington- 
Finkclstein coordinates [12j | . For our approximate wave- 
forms decribed in § [III m(u) is M constant outside, it 
is quickly changing within a short range < u < Ait-tot 
across, and it is M/ = Mq — Am to t inside the radiation 
shell. Here, Am to t and Au to t are set by the simulated 
waveforms § [III 

For the metric given by Eq. ([5]) it is straightforward 
to derive the convergence of null-geodesics describing the 
world lines of radiation shells using Raychaudhuri's equa- 
tion or the equation of geodesic deviation 68] . In either 
way, we find that radially outgoing shells of radiation sim- 
ply follow the world lines u(r) = constant, implying that 
the Am coordinate difference between the shells does not 
change during the propagation. The physical contrac- 
tion of radiation shells can be examined using the proper 



time measure between shells and the observed luminosity 
profile. 



A. Proper time duration 

First, we estimate the proper-time duration of the GW 
signal along the world-line of a hypothetical observer 
crossing the radiation shell. For simplicity, we restrict 
to observers at a fixed spacial coordinate (r, 9, <p). Then 
we have dr = d9 = dfi = and Eq. §5§ gives 



dr 2 



-ds 2 = 1 - 



2m(u) 



du 2 



(6) 



leading to dr = y/l — 2m(u)/rdu. Therefore, du can 
be interpreted as the infinitesimal proper time difference 
between two radiation shells at fixed radius approaching 
infinity. Thus, we adopt the notation du = dr^, and 
similarly Au = At,^ for integrated quantities. Finally, 
let us define m t = m(ui) for i € {1, 2} and Am = \m,2 — 
m\\. We set (mi, m-2, Am, to t) = (Mo,Mf,eMo) when 
referring to the total GW signal duration. 

Integrating between two arbitrary shells of radiation 
U2 and ui gives 



At = dr = 



1 - 



2m(u) 



1/2 



du. 



(7) 



After expanding m(u) in a Taylor-series, the integrand 
becomes 



1 - 



2mi 2u(m'} u 2 (m") 



Substituting in Eq. ((7]), and using (m') 
and (m") = to first order, we get 



1/2 

• (8) 
-Am/A-Too 



Ar(mi,m 2 ,r) 



rAr^ 



1 - 



2m 



3/2 



(9) 



3|m2 — mi 

Setting mi = M and expanding in terms of Am, we get 
At 



1 



2M 



1 



1 Am 



2 r - 2M Q 



0{Am 2 ) 



( 10 ) 

The leading order term can be indentified as the gravi- 
tational redshift for constant mass Mq, the second term 
describes the correction due to the radiation mass. If we 
expand also in terms of powers of l/r, the relative change 
in the proper time duration of the signal becomes 



At - A? 

AToo 



= - M, 



Am\ 1 

— )- + 0(r- 2 ,Am 2 ). 11 
2 I r 



Equations (IMTTj) describe the self-gravitational distor- 
tion between radiation shells in terms of proper time 
along world lines of r — const. One can notice that to 
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leading order this is simply the gravitational redshift for 
the average enclosed mass between the shells. However, 
since Am changes along the wave packet non-trivally for 
a fixed radius as a function of time, the modification of 
the profile is generally not self-similar, leading to the dis- 
tortion of the luminosity profile as a function of radius. 

To make this point clearer, we correct for the average 
distortion of the signal and calculate the residual distor- 
tion of the signal. Let us define 



At' = (1 + z)At - A? 



(12) 



where (1 + z) is a time-independent constant representing 
the "average gravitational redshift" at a given r, given by 



l + z = 



1 



1 - 



2{m) 



(13) 



(m) = 
order 



(mi + m,2)/2 is the average mass, and to leading 



(m) _ M 
r r 



Am 
~2r r ' 



(14) 



Now let us take an arbitrary radiation shell enclos- 
ing mass Am relative to the shell enclosing mass Ma — 
0.5Am to tj i-e. we set mi = Mo — 0.5Am to t and m^ — 
mi — Am in eq. ©■ After correcting for the average grav- 
itational redshift using Eqs. (fl~2]) and (fl3|) . the residual 
relative distortion to leading order is 



At' Am 9 . 



(15) 



Figure [2] shows the residual distortion using the exact 
formula (Eq. (fT2|) , thick lines) and the leading order con- 
tribution (given by Eq. (fT5")) . dotted lines). At the typ- 
ical radius used by numerical simulations for waveform 
extraction, r — 50Mq, the primary bulk waveform distor- 
tion changes the signal duration by M /r ~ 2%, and the 
secondary relative waveform distorsion between the front 
and the back of the signal is eMo/(2r) ~ 7 x 10~ 4 for a 
typical BH inspiral-merger with high spins e = 7% (see 
§ |TT| . The figure also shows that the higher order effects 
beyond 1/r lead to an uncertainty of order 10 _3 -10~ 4 for 
r = (30-50)M . 



B. Luminosity Profile 

Since Eqs. © and (fTTj) are applicable to two arbitrary 
shells of radiation, we can use them to compute the evo- 
lution of an arbitrary initial radiation profile, whereas the 
luminosity is simply L = Am/ At, the total mass-energy 
crossing a sphere at radius r within proper time At. 

The profile at infinity is given by m(u) in radiation 
coordinates, or At^to), the proper time a shell enclos- 
ing mass m arrives at r = R where R — > oo, relative to 
the outermost shell of radiation [84[ . Here A t^ (m) can 
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FIG. 2: The residual self-gravitational distortion of shells af- 
ter correcting for the bulk gravitational redshift. The thick 
curves show the change in the proper time duration of the 
signal at radial distance r from the source between shells en- 
closing 7% (top) or 3% (bottom) of the total mass, dotted 
lines correspond to the leading order term given by Eq. Q5) . 
The vertical lines highlight the typical radii used in numerical 
simulations for GW extraction. 



be any monotonically decreasing function, for which the 
luminosity at R in Eq. ([3]) is 



Loo(t) = - 



d Atqo (to) 



dm 



(16) 



where the minus sign originates from our definition of m: 
the shell labeled by the largest value of m arrives the ear- 
liest. The luminosity profile can also be obtained as the 
function of time, Loo(t), using the relationship At(to). 
Conversely, for given Loo(t), we can compute Atoo(to) 
using Eq. (fTo) . The luminosity at some other distance r 
can be obtained similarly if given Ar(m, r), the arrival 
time of mass m relative to the outermost shell at distance 
r. This function is given by Eq. ©, substituting the 
waveform AToo(m) for Atoo, and (toi,TO2) = (Mo, to). 
The luminosity profile at r is then 



L r (m) = — 



9AT(m, r) 



dm 



2m 



-1/2 

1 - — I Loo(m) 
r J 

(17) 



Therefore, the modification of the profile in Bondi ra- 
diative coordinates (m, r) the profile is distorted self- 
similarly. However, in terms of the observer proper time 
variable, Ar(m) = f™ L r (m)~ 1 dm, L r (r), the modifica- 
tion to the profile will not be self-similar: 



L r (r) 



2f Q T L r (T>)dT r 



-1/2 



L r (T / )dr / 
(18) 
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Equation (|18p relates the luminosity profile as a function 
of proper time at radius r, L r (r), to the profile at infinity, 
ioo(r). Comparing Eqs. (fT7|) and (fl~8|) shows the advan- 
tage of Bondi type radiative coordiantes as opposed to 
proper time. 
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FIG. 3: Our fits to the GW luminosity profile for binary BH 
inspiral merger simulations as a function of observer proper 
time t and the evolution of the profile at various distances, r. 
The profile parameters are given in §|II]and e = 7%. Top: The 
absolute profile (in units of c 5 /G) is shown for two extremes, 
a nearby distance (r = 12.5Mo) and far-away distance (r = 
10 3 Mo). Bottom: The difference between the GW luminosity 
profiles at infinity (i.e. r — 10 3 Mo) and three cases of smaller 
r, in units of peak luminosity at infinity. The peaks of the 
profiles are set to t = 0. The main effect responsible for the 
differences seen in this figure is the bulk gravitational redshift. 

Figure [3] plots L r (r) for our fit to the luminosity profile 
at infinity of merging binary BHs Loo(t) with e = 7% (see 
§|TT|. The top panel shows the absolute profile while the 



bottom panel shows the difference between the profile at 
some radius r and the profile at infinity, such that the 
peak of the profiles are at r = 0. The bottom panel 
is useful to visualize the characteristic evolution of the 
profile. 
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FIG. 4: The residual self-gravitational distortion to the lumi- 
nosity profiles after accounting for the avarage gravitational 
redshift, z = Mo/r (top) or (Mo — 0.5Am to t)/r (bottom), re- 
spectively. Other parameters are the same as in Fig. 3. 



In § IIII A| we have identified the two main effects re- 
sponsible for the convergence rate of the waveform to 
be the gravitational redshift corresponding to the aver- 
age mass and the self-gravitational effect. Indeed, the 
differences visible in Figure [3] are primarily due to the 
former. However, correcting for only the average gravi- 
tational redshift at each radius r leaves a nonnegligible 
systematic error with respect to the true signal. To see 
this, we substitute Ar 00 /(1 + z) given by Eq. (fT3"|) into 
Eqs. (|16[) and (|17p. and refer to the corresponding lumi- 
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nosity as the average gravitational redshifted luminosity 
profile, L z r (t) . After subtracting from the true profile for 
each t, the residual luminosity distortion is 

L' r (r) = L r (r) - L*(t), (19) 

where we set the reference time again to r = for the 
peak of the luminosity profiles. 

Figure 3] shows the residual self- gravitational distortion 
L' r (r) for various radii in units of the peak luminosity at 
infinity, Xoo.max- Naturally, the definition of the "av- 
erage gravitational redshift," z, used for defining L' t {t) 
makes a difference in the result. The top panel uses only 
the initial binary mass (m) = Mg/r in Eq. (fT4]) totally 
neglecting the gravity of the radiation shell, while the 
bottom panel has (m) = (Mo — 0.5A?7Jtot)A, i.e. the 
redshift is chosen to account also for the average gravity 
of the radiation shell. In the later case we find a much 
quicker convergence for the waveform peak at increasing 
radii, but the former choice is more suitable for the early 
parts of the waveform corresponding to the late inspiral 
waveform. A comparison of Fig. [3] and 2] shows that the 
self-gravitational distortion is roughly an order of magni- 
tude smaller than the effect of the average gravitational 
redshift. 



C. Self-Gravitational Coordinate Effects 

In the previous sections we have derived the time du- 
ration and the luminosity profile of the radiation shell as 
it propagates radially outward from the source. We have 
assumed that the GW profile at each fixed arial radius r 
is parameterized by the proper time r of a hypothetical 
observer fixed at that radius, in particular the luminosity 
L r (r) was the total mass-energy crossing a sphere at ra- 
dius r within infinitesimal proper time dr. Therefore, the 
adopted time-coordinate r corresponds to a synchronous 
gauge at each radius. Since physical observables depend 
precisely on proper measures, these coordinates allow a 
simple interpretation of the convergence characteristics 
of the GW profile at large radii. 

Other choices of coordinates would have introduced 
additional artificial distortion effects making the conver- 
gence characteristics of the waveforms much different. 
Consider for instance the "natural" coordinate system 
(t, r, 6, 4>) in the spherically symmetric case that is chosen 
to be Schwarzschild both before and after the GW burst 
has arrived with masses Mq and Mf = Mo — Am to t, re- 
spectively, and which changes smoothly in between these 
regions. An example of such a coordinate system can be 
derived from the Vaidya metric Eq. ([5]) with the implicit 
transformation t = u + r + 2m(u) ln[r — 2m(u)] where 
m(u) describes the mass function interior to u (which is 
constant along the outgoing radiation world lines). In- 
deed, everywhere in the spacetime where dm/du = 0, 
these coordinates yield a Schwarzschild metric, and the 
Vaidya metric in radiation coordinates ^ is then sim- 
ply the Schwarzschild solution in Eddington-Finkelstein 



coordinates [12| in these regions. This map covers all 
relevant parts of the spacetime including the GW zone. 
The world-lines of radiation shells can be shown to follow 



where the second term is called the Shapiro-time delay 
\6l for a particle crawling out of the gravitational poten- 
tial of mass m u interior to it. After integration, we find 
that to first order the temporal separation of two radi- 
ation shells enclosing mass Am at radius r evolves to 
leading order as At(r) — At(Ro) — 2Amln(r / Rq), where 
Ro is an arbitrary initial radius. In these coordinates, 
the signal duration contracts uniformly in exponential 
distance intervals. Even though the metric is asymptoti- 
cally Minkowski (where t approaches r for r 3> Mo), the 
resulting profile evolution is fundamentally different from 
Ar(r) given by Eq. (fTTj) ! 

The appearance of the logarithmic radial dependence 
of the waveform was first realized by Fock . This ef- 
fect is specific to the harmonic coordinates and can be 
avoided if changing to Bondi type radiative coordinates 
[70l . PrO ] . Blanchet & Schafer [38| have shown that a sim- 
ilar logarithmic dependence of the GW tail leads to a 
tail-induced amplitude and phase shift (typically of order 
1CP 7 ) for stationary sources. In contrast, the logarith- 
mic radial dependence of the wave contraction for merger 
waveforms can be significant for GW merger simulations. 
Between r = (20-40)Mo, the waveform contracts in At 
by a fraction of 5 x 10 -3 , which is just of the order of the 
current wave extraction precision 0, [26|, [2tJ • This appar- 
ent logarithmic contraction effects can be avoided if one 
changes to the proper time variable as we have done in 
the previous sections. The remaining Ar(r) evolution is 
however a physical effect. 

Both the logarithmic At(r) contraction and the phys- 
ical Ar(r) evolution (in particular the contribution de- 
noted by Ar'(r) above) are consequences of higher order 
radiation effects in the Einstein equations beyond the 
scope of first order methods such as the Regge-Wheeler- 
Zerilli-Moncrief perturbation method used for extrapo- 
lating the numerical waveforms to infinity. Therefore 
these effects cause the extrapolated waveforms to be 
different when extracting GWs from numerical simula- 
tions at various radii by standard methods using no self- 
gravitational interaction. For a related recent analysis 
see Ref. [72j]. 

IV. DISCUSSION 

A. Summary 

We considered the self-gravitational effect of gravi- 
tational radiation on the propagation of GWs from a 
compact source. We adopted simple approximations for 
the geometry of the radiation, by considering spherical 
symmetry on scales comparable to the radial width of 
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the radiation packet. This approximation appears ade- 
quate for the quadrupolar (I = 2, m = 2) radiation pat- 
tern around binary BH sources in numerical simulations 
[25l l26l . I27I ]. Nevertheless, we use the Appendix to ex- 
amine the maximal effects of anisotropy in the opposite 
(exaggerated) extreme, when the outgoing radiation is 
concentrated into a compact region. We find that irre- 
spective of the level of anisotropy, the gravitational radi- 
ation is distorted under the influence of its own gravity 
as it propagates. Contrary to the standard gravitational 
redshift, which is a uniform shift of the waveform, the 
self-gravitational effect depends on the intensity and is 
predominant only for the most intensive bursts of radi- 
ation causing a non-uniform distortion of the waveform. 
The self-gravitational distortion depends on distance to 
leading order as Am/r, and is therefore relevant on scales 
r sg /Mo <e/6 where e is the radiation efficiency and 5 is 
the desired calculation accuracy. For BH binary merg- 
ers simulations e ~ 7% and 5 ~ 10~ 5 implying that 
r sg < 7 x 10 3 Afo- If the GWs are extracted within this 
region, the self-gravitational distortion should be taken 
into account. 



B. Testing the Effect with Numerical Simulations 

Numerical simulations based on the full set of Einstein 
equations for binary BH inspirals have not yet reported 
evidence for the waveform distortion effect considered 
here although they have shown that the waveforms do 
not converge within a fractional accuracy of S ~ 10~ 3 
[S HI, H3] • This is because the simulations are restricted 
to a limited volume, typically of radii ~ (80-850)Af, 
while the extracted waveforms are typically compared be- 
tween r ~ 20-50M. For a radiation mass Am/Mo ~ 7%, 
the primary effect is a shift of the waveform due to a 
logarithmic Shapiro time delay of the remnant, a uni- 
form gravitational redshift, and the self-gravitational ef- 
fect. We have shown that the logarithmic Shapiro de- 
lay does not show up if using proper measures to de- 
scribe the waveform, and the uniform gravitational red- 
shift is accounted for in the linear wave propagation mod- 
els. However, the residual self-gravitational effect in the 
GW luminosity has a characteristic profile that has to be 
subtracted when extrapolating the extracted waveform. 
The peak of the effective luminosity distortion reaches 
2 x 1(T 3 and 5 x 10~ 4 at r = 30 and 50M , respectively. 

Present-day numerical relativity simulations should al- 
ready be capable of directly measuring the relevance of 
our effect by artificially amplifying the gravitational radi- 
ation found at the extraction radius r ~ 20 M, and start- 
ing the simulation with these amplified initial conditions. 
For example, for a total GW energy Am/Mo = 30%, the 
effective luminosity distortion between r = 20M between 
20-50M is several percent, which is well within simula- 
tion and extraction errors. For consistency, the simu- 
lation should confirm that the total energy content of 
the radiation does not change. We also expect the initial 



ringdown frequency (corresponding to the most energetic 
shell) to be smaller than the final ringdown frequency. 

In order to avoid errors caused by the self-gravitational 
distortion effect up to the desired numerical precision S ~ 
10~ 5 , the waveform extraction radius should be chosen 
to be r sg > 7 x 10 3 Mo- Alternatively, if the waveforms 
are extracted at smaller radii, the waveforms should be 
converted to Bondi type radiative coordinates and then 
extrapolated with the scaling 1/r. 



C. Observational Implications 

The self-gravitational waveform distortion is important 
for future observations of BH binary mergers. 

1. The waveform distortion is expected to be resolv- 
able for the LISA instrument with respect to sim- 
ulated waveforms for total BH masses of (10 4 - 
1O 9 )M0. The total signal to noise ratio of merger 
waveforms is 10 4 for LISA observing z c ~ 1 |. 
The distortion effect modifies the waveform ampli- 
tude and frequency by ~ (10 _3 -10 -4 ) for numerical 
waveforms extracted between r = (30-80)Mq. 

2. The distortion involves a systematic modification 
of the waveform which needs to be accounted for 
in order to interpret observed merger waveforms 
and improve the estimation uncertainty of physical 
parameters beyond the uncertainty of the preced- 
ing inspiral signal. The signal to noise ratio of the 
final BH merger waveform is an order of magni- 
tude larger than for the inspiral, implying that the 
merger waveform has a potential to greatly reduce 
parameter estimation errors. Note that the relative 
accuracy using only the inspiral sig nal with LISA 
is expected to be 10~ 3 -10~ 5 [73],[74[ for estimating 
the component masses, which is smaller than the 
distortion effect. 

3. This effect is different from the uncertainties caused 
by gravitational lensing (7f|, in that it is only an 
issue concerning the convergence properties of nu- 
merical simulations. Lensing causes an error of 
several percent on the inferred luminosity distance 
(due to the unresolved matter along the line of 
sight), and lensing errors increase with the source 
distance. In contrast, the self-gravitational effect is 
of order 0.1-0.01 percent for numerical simulations 
if the waveforms are extracted at 30M-80M and 
dies off quickly as 1/r. The wave distortion effect 
is of order 10 -20 relative to the waveform ampli- 
tude for typical astrophysical scales. Therefore, the 
wave-contraction effect does not provide any addi- 
tional physical parameters for observations. 
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APPENDIX A: ANISOTROPY 

Our analysis considered only perfectly spherically- 
symmetric configurations. The approach was motivated 
by the quasi-spherical radiation patterns found around 
binary BH sources in numerical simulations [25|, [26|, [2tJ • 
In this Appendix, we would like to examine the sensitiv- 
ity of our basic results to deviations from sphericity. To 
gauge whether there is any such sensitivity, we analyse 
the most extreme case in which the outgoing radiation is 
concentrated into a highly compact region. 

But first let us define more precisely what we assumed 
so far. The derivation presented in §|TTT] requires that the 
radiation field is "initially locally spherically symmetric" , 
so that it is initially described by the Schwarzschild met- 
ric locally within some narrow solid angle AO < Ar /r be- 
fore the GW arrives, where Ar ~ cAt is the radial width 
of the wave-packet along its propagation direction. But 
since the radiation propagating through this solid angle 
is in no causal contact with the radiation field expanding 
towards other directions, it cannot distinguish the ac- 
tual spacetime from a spherically symmetric one. Note 
that the outermost shells of radiation expanding along 
different directions are always causally disconnected by 
definition, and the interior shells of radiation can only be 
affected by the outer shells within A9. Since we examine 
the distortion effect on large distances compared to the 
width of the burst r 3> Ar, spherical symmetry must 
only be required within a very narrow angle. This simple 
set of considerations implies that if high-order multipoles 
have a vanishing contribution at large distances, the re- 
sults derived in § IIIII are applicable very generally for 
short bursts of radiation, Ar <C r. Indeed, numerical 
simulations confirmed that the dominant contribution to 
the wave amplitude is given by the (I = 2, to = 2) multi- 
pole and higher order terms are suppressed by more than 
a factor of magnitude (see references in § [IT| . In the re- 
mainder of this Appendix we demonstrate the validity of 
this simple conclusion through explicit calculations. 

We consider three variations on a toy model to estimate 
the effect of anisotropy. We start with the simplest model 
and refine this model by adding more details and com- 
plexity in the successive models. In each case, we discuss 
general implications for the model under consideration. 
In all models we consider the extreme opposite regime to 
spherical symmetry, namely that the radiation is maxi- 
mally clumped into two outgoing BHs L (leading) and T 
(trailing) of masses mi and tot, representing the leading 
and trailing edges of the radiation, respectively. We as- 



sume that L and T are moving in the same direction on 
light-like world lines, so that T lies in the causal past of 
L, but L is outside the causal past of T throughout their 
propagation. We assume that there is also a remnant 
Schwarzschild BH R with mass rriR. The instantaneous 
radial position coordinate of R, T, and L at time t are 
0, rr(t), and rz, (t). We are interested in obtaining the 
world lines of BHs T and L to see how the coordinate 
separation Ar(t) = 7*1, (t) — rr(t) decreases with time as 
compared to the spherically symmetric result. In our 
first model we neglect the remnant R (setting = 0), 
assume that L moves with constant velocity Vl in free 
space, and calculate the trajectory of T in the spacetime 
created by L. Subsequently, we will generalize L to move 
on a more general world line with a slowly changing ve- 
locity Ui(t). Finally, we can turn on the remnant R in 
addition to L, and include the retardation effect when 
calculating the relative motion of T. 

We note that the spacetime of BHs moving at the 
speed of light have been calculated previously in Ref. , 
which found that BHs moving in the same direction do 
not interact. However, Ref. [46[ assumed that the BHs 
move in free space and consequently adopted v = 1 for 
their velocity. In contrast, the BHs T and L travel on 
null-geodesics in the perturbed spacetime which is ini- 
tially the Schwarzschild spacetime. This difference gives 
rise to a non-trivial interaction between the BHs T and 
L. 

We compare our results to the spherical case, using the 
(t, r) coordinate system defined by Scwarzschild coordi- 
nates before and after the radiation shells as described in 

§Ena 



1. No remnant rriR — 0, constant vl velocity 

We start by assuming that L is a BH with constant 
velocity vl < 1 in free space, and wish to calculate the 
world line of T in this background. Here, we assume that 
no remnant is present, and that L and T move along the 
same spatial direction, which we denote by x. Thus it is 
sufficient to restrict our attention to the two dimensions 
(i, x) of the spacetime. 

Let us start by deriving the metric. In the coordi- 
nate system (t',x') comoving with L, the metric is the 
Schwarzschild metric ds 2 = -(l-(j)')dt' 2 + (l-(j)')- 1 dx' 2 , 
where </)' = 2m' L /\x'\. Here x' = corresponds to the BH 
L for all t', and m' L = m L /j is the rest mass of L, where 
tol is the energy carried by L in the original (t, x) co- 
ordinates and 7 = 1/a/I — v\ is the Lorentz factor. To 
derive the metric in the (t, x) coordinate system, we ap- 
ply the diffeomorphism (t,x) = j(t' + vr,x' ,vi,t' + %'), i.e. 
a global Lorentz transformation, 



ds 7 



2v L cj ) '(2 - ft) 

(1-t,£)(W) 



-dt l 



1 - v 2 L (l - cj>') 2 

(l-iiXW) 



dtdx 



dx 2 



(Al) 
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which can be rearranged as 



ds 2 



[v L + (1 - <j>')]dt - [1 + v L (l - </>')]dx 

V(i-t£)(i -0') 

Jv L - (1 - </>')]dt - [1 - V L {1 - </>')]dx 



(A2) 



Here is to be expressed as the function of the new 
coordinates (t,x), i.e. 0' = 7~ 2 where </> = 2mz,/|Ax|, 
Ax — xl ~ x, and x^ = VLt is the instantaneous position 
of the singularity. 

The xt it) null-geodesics describing the world line of 
T can be obtained by setting ds 2 = 0. Equation (|A2|) 
shows that there are two solutions 



[vl + (1 - 7 _2 0)]dt - [1 + v L (l - ~f- 2 0)}dx T = 0, 



(A3) 

[v L - (1 - 7~ 2 0)]d* - [1 - v L (l - 7 ~ 2 0)]dx T = 0. 

(A4) 

These differential equations can also be obtained more 
simply by finding the null-geodesics in the comoving co- 
ordinates (t',x') first, and changing to the (t,x) coordi- 
nates only in the resulting equation. The null geodesies in 
the comoving coordinates are simply dx' T /dt' = ±|1 — <f>'\ 
(see Eq. [20]), and the Lorentz boost coordinate trans- 
formation of this differential equation leads instantly to 
(IA3HA4|) . Therefore, the two solutions (|A3flA4|) describe 
the null geodesies approaching or receding the moving 
BH, respectively. 

We would like to find the solution for the T test particle 
approaching the source L from behind, namely Eq. (|A3[) 
for an initial condition xt < Xl- This first-order dif- 
ferential equation can be solved analytically by a linear 
substitution. The coordinate velocity vt = dxT(t)/dt 
monotonously decreases from 1 to vl as the event horizon 
at XLhor(t) = VLt — 2j~ 2 itil is approached. In particular 
if Vl = 1) i.e. the source L has the speed of light in the 
free-space background, then the trailing test particle T 
will not be delayed at all, VT{t) = 1 for all t. However, 
if Vl "C 1 then T is considerably affected by the Shapiro 
delay near the horizon of L. 

In concluding the description of this model, let us sum- 
marize how the clumpy case compares to the spherically 
symmetric case of expanding radiation shells. First re- 
call that in the spherically symmetric case, L has no ef- 
fect on T throughout the dynamics regardless of vl or 
Ar. In the clumpy case, the gravity of L delays the mo- 
tion of T. The magnitude of this delay is significant only 
if both of two conditions are violated: (a) vl ~ 1 and 
(b) Ax = xl — xt ^> 2-f~ 2 rriL- What are the "typi- 
cal numbers" for these quantities? Eq. ([2"0")) implies that 
Vl = 1 — 2Mq/x (which is also true in the clumpy case, 
see § IA3|) . implying that 7~ 2 ~ 4M /r, and for binary 
mergers Ax ~ 12ttM , m L < M - M f < 0.06Af , we 
find that the two cases are equivalent to (a) x ^> 2Mq 
and (b) x > 0.2Am > 0.01M . Quite clearly, these 
conditions will not be violated for distances outside the 



dynamical regime of strong gravity e.g. x > Rq — 30Mo. 
Thus, we expect only very minor modifications relative to 
the spherically symmetric case, even in the most clumpy 
case. For a quantitative estimate we need to integrate 
these modifications over the relevant distances which we 
describe next. 



2. No remnant, slowly changing «z,(t) 

Next we consider a slowly changing source velocity 
fi(i) for the BH L, continue to neglect a remnant R, 
and calculate the motion of T in this spacetime. Since L 
is assumed to move on a light-like world-line, it is not ef- 
fected by R and only responds to the background created 
prior to the production of the bursts. Thus we assume 
L moves on the null-geodesic of the background as de- 
scribed by ([20)1 with vl = 1 — 27tit/xl- 

If vl is slowly changing, we can consider vl to be con- 
stant during short time intervals with infinitesimal jumps 
on their boundaries. We can then find the correspond- 
ing world line segments of T by solving the differential 
equation (|A2p and matching the boundary conditions of 
the successive segments by requiring continuity. In the 
limit that the length of the constant time intervals ap- 
proaches zero, the world line xx{t) at every instant is 
given by Eq. (|A2|) with vl now denoting the instanta- 
neous velocity, and <fi' referring to the instantaneous value 



of the potential: 



7" 



Aa; = xl — x and xl — J VL{t)dt. Thus, 



where <f> — 2niL/\A.x\ with 



dxT 

~df " l + v L (l 



v L + l- 7~ 



T 



2m L (l-v L y 



Ax 



2rn L v L {l - Vl) ' 
(A5) 

Substituting Vl = 1 — ^itit/xl, the distance between the 
clumps of radiation satisfies 



dAx 2rriT 
~dT = 



ArjiLrriT 



XL 



xlAx — 4mimT (l — 2mr 




Note that xi,(t) can be used to express the width of the 
packet Ax as a function of xl- To simplify the result, 
we use qL = m/,/my, set the units to the Schwarzschild 
radius 2niT — 1, and express (|A6p in terms of the loga- 
rithmic distance variable, y = ln(xL — 1). Then 



dAx 
dy 



= -1 



qL 



x L Ax - q L (l 



which to first order in 1/xl becomes 
dAx _ qL 



dy 



-1 



x^Ax 



(A7) 



(A8) 



Equation (|A8|) shows that to leading order, the wave- 
packet packet contracts linearly in terms of the logarith- 
mic distance variable x. In the spherically-symmetric 
case, the inner shell T is not influenced by L and so 



12 



dxr/dt = 1 instead of Eq. (|A5[) . leading to dAx/dy = 
— 1. Therefore the distortion of wave-packets in the max- 
imally clumpy case is the same as in the spherically sym- 
metric case to leading order. The difference arises in the 
next order given by the second term in (|A8[) describing 
how the gravity of L Shapiro-delays the motion of T. 



from which 



This is typically of order (4 x 6%)/(127r) x x i 



io- 



for xl ~ 30Afo and gets exponentially smaller for expo- 
nentially larger distances. 



3. Remnant included, changing vl 

The previous model assumed no remnant (i.e. m# = 
0), and postulated that T propagated in the spacetime of 
L by assuming that the spacetime at T was the spacetime 
of L at the same instant. Here we consider a nonzero 
mn and account for the retardation of the effect of L as 
percieved by T. We assume a slowly changing velocity 
and that the gravitational perturbations are sufficiently 
small to allow simple superposition to leading order. We 
follow a simplified approach with the following essential 
assumptions: 

1. The initial condition is a spherically symmetric 
Schwarzschild spacetime centered at x = 0. 

2. L moves on a null- geodesic Xz,(i) of the initial back- 
ground metric of R and T, i.e. in a Schwarzschild 
metric centered at x = for all t and for a mass 
mji + tot- The world- line of L follows ([2H|) accord- 
ingly. 

3. T moves on a null-geodesic Xr(t) of the background 
metric of R and L, which we assume to be a sim- 
ple superposition gij = r]ij + Sg^j + Sgfj' ret . Here 
ddij = 9ij - Vij f° r a given metric, gfj , where 

is the Minkowski metric, gf^ is the metric of the 
remnant i.e. the Schwarzschild metric centered at 
x = for all t with mass rriR. The metric ghjs the 
stationary boosted Schwarzschild metric (|A1HA2[) 
with mass rriL, velocity vl, centered at xl- The 
label rct stands for retardation, which we describe 
next separately. 

4. We account for retardation by setting At lct = t — 
t re t to be the light-travel time from L to T. For this 
we compute the inward propagating null-geodesics 
from L to T [i.e. between positions (i ro t, ^i(iret)) 
and (t, rT(t))] , based on the initial background met- 
ric of R and T (i.e. neglecting the gravity of L). 

To find the retardation time, we note that the inward 
propagating null geodesies satisfies dx/dt = —[l — 2(mR+ 
tot) / x] . Since this is exactly the time-reversed world line 
of L, we get XT(t) = xr.it — 2Ai ret ). Integrating dt/dx 
for the world line of L between Xi(t — 2Ai rct ) and XL(t), 



2At let = 



x L (t) 



XL 



x T (t) x l - 2to^ - 2m T 



dx L , (A9) 



Ai rc 



XL - X T 



+ (ttir + m T ) In 



xl — 2mft — 2tot 



xt — 2m a — 2tot 

(A10) 

The distance where T percieves L is x^it) — XL(t — 



At re t) and the separation is Ax rc t 



xt- Substi- 



tuting xt as xt — xl ~ Ax and xt = x T £ — Ax ret into 
(|A10() . we can find Ax Ic t for given Ax and xl- To first 
nonvanishing order in 1/xl, 



Ax rR t — 



[rn R 



tut) Ax 2 



4x1 



(All) 



The leading order term corresponds to the propagation at 
the speed of light in free space, vl = Vt = 1- Note that 
the correction is proportional to x^ 2 , which is extremely 
small for the physical cases beyond the strong field zone. 
Finally, we define the retarded position and velocity 
xf* = x L - Ax Tet , and v 1 ^ = 1 - 2(to^ + m T )/x T L st , 
which can be written in terms of Xl and Aa; using Eq. 

HMD) . 

The metric contribution Sg^' let of L at T at time t, is 
the boosted Schwarzschild metric (|AllA2p with instan- 
taneous velocity v L ct , a singularity at x™*, and distance 
Ax rct . 

Now we can redo the derivation presented in § IA 21 to 
find the motion of T, using the modified spacetime g^ 
given above. Again we find two solutions for ds 2 = rep- 
resenting the ingoing and outgoing radiation. Expanding 
the outgoing solution in a series in x^ 1 , we find 



dxT 
~df 



clumpy 



2te 

XL 



QrAx 



2q L (l + q 2 R ) 



Ax 



1 



(A12) 

where qi = to,/(to_r + tot) for i € {R, T, L} and distance 
units are chosen to be the Schwarzschild radius 2(mn + 
tut) = 1- In order to get the instantaneous shell width 
Ax as a function of logarithmic distance y, we can redo 
the manipulations of l|A6t4A8p for the result (|A12|> . To 
first order in 1/xl, 



dAx 



dy 



-q T 



clumpy 



qnAx 



2q L (l + q 2 R ) 



Aj 



1 



XL 

(A13) 

In the limit of no remnant qn = 0, we almost recover 
the solution derived previously in Eq. (|A8|) . There is a 
factor 2 difference, which is the direct consequence of the 
retardation of the percieved distance Ax, which had been 
neglected in Eq. (|A"8|) . 

Equation (j A13|) should be contrasted to the expansion 
of two spherically symmetric shells tot and in the 
presence of a remnant m R . We may expand the corre- 
sponding spherical solution of § IIII CI in a series in 1/xl 
to first order: 



dAa 



dy 



= -qr 



spherical 



XL 



(A14) 
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The first two terms in Eqs. (|A13[) and (|A14[) are iden- 
tical. The correction describing the "Shapiro delay" of 
contraction in the clumpy case due to the gravity of L is 
2<7l(1 + qj t )/(xL^x). Substituting typical physical val- 
ues qn = 97%, qx = 3%, qn — 3%, and Ax Q ~ 6tt ~ xlo 
in units of Schwarzschild radii, we see that the correc- 



tion is of order 10 -4 initially, and becomes exponentially 
smaller at exponentially larger distances. In summary, 
even in the most extreme case of dumpiness, the radia- 
tion packet propagates to very high precision according 
to the spherically symmetric description. 
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r — > oo is the Minkowski metric, neglecting cosmologi- 
cal effects. This is justified by the orders of magnitudes 
or radii where the radiation shell distortion is relevant 
r < 10 3 A/o (see Fig. [2j is much less than cosmological 
scales r > 10 18 M . 



